Band-gap structure and singular character of bounded 
one-dimensional multibarrier potential 



D. Bar 



February 1, 2008 



Abstract 

The bounded one dimensional multibarrier potential shows signs of chaos, phase transition and 
a transmission probability of unity for certain values of its total length L and the ratio c of total 
interval to total width. Like the infinite Kronig-Penney system, which is arranged along the 
whole spatial region, the bounded multibarrier potential has a band-gap structure in its energy 
spectrum. But unlike the Kronig-Penney system, in which the gaps disappear for large energies, 
these gaps do not disappear for certain values of L and c. The energy is discontinuous even in 
parts of the spectrum with no gaps at all. These results imply that the energy spectrum of the 
bounded multibarrier system is singular. 
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1 Introduction 



It is known that the Kronig-Penney multibarrier system ^ [2] which is arranged along the whole 
spatial region is characterized by an energy spectrum which is composed of continuous bands 
separated by forbidden gaps |Sj . The larger are the values of the energies the shorter these 
gaps become until they entirely disappear for large enough energies It is also known 

that the gaps in the energy spectrum of the infinite Kronig-Penney system are related to the 
eigenvalues A± of the characteristic equation which is associated with the two dimensional transfer 
matrix that relates the transmission and reflection coefficients of any two neighbouring barriers. 
The limits (linv^-i-oo A!^) are taken when the related potential barriers are separated from each 
other by a very large number n of other similar barriers. Thus, if either of these limits becomes very 
large, in which case the resulting wave function does not remain finite, the corresponding energies 
are disallowed [T] and these constitute the gaps of the energy spectrum. 
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The former problems which are associated with the infinite Kronig-Penney system do not arise in 
the bounded multibarrier system discussed here. This is because although the number of barriers is 
very large, as for the Kronig-Penney system, nevertheless the finite extent of the system enables one 
to analytically express in explicit form [3J Ej the total transfer matrix that relates the two barriers 
at the two extreme sides of the system. Thus, the characteristic equation and the corresponding 
eigenvalues of this overall transfer matrix have also been found [2j in exact closed forms. From 
the latter one may construct the appropriate energy spectrum for both cases of E > V and E < V 
where E is the energy and V the constant height of all the barriers. 

An important parameter which will be constantly referred to in this work is related to the mentioned 
eigenvalues X± and is the variable over which the energy E is plotted (see, for example, Figure 6.6 
in pQ and the seven figures here). This parameter is denoted n here and is defined by the first of 
the following Eqs (see also Eqs (|5|l-(|lU|)). We show in the following that there are no allowed 
energies that correspond to certain values of the total length L of the system and the ratio c of its 
total interval to total width. Moreover, in contrast to the Kronig-Penney system in which the gaps 
in the energy spectrum disappear ; l a for either large values of the energy E or (and) large values of 
the parameter n [I] the case here is different. This is because, as will be shown in Sections 3-4, the 
gaps in the energy spectrum of the bounded one-dimensional multibarrier system do not disappear 
for certain values of L and c even for large E or (and) large k. We also show that the energy E 
depends upon the total length L and the ratio c in such a manner that its form, as function of k, 
is quasi-periodic for small c and large (or intermediate) values of L and is proportional to k for all 
c > 3. It will also be shown that the energy spectrum is discontinuous at many of its points even 
at the parts at which it has no gaps at all. 

In Section 2 we introduce the terminology and terms [H that describe the bounded multibarrier 
potential. In Section 3 we discuss the energy spectrum for the E>V case and find the values of the 
parameter k for which no corresponding energies are found. That is, we find the band-gap structure 
of it. We also show the remarked discontinuity in E as function of k. In Section 4 we repeat the 
whole process for the E < V case and summarize in Table 1 the allowed energies for both cases of 
E>V and E < V. In Section 5 we show that the occurence of stable gaps in its energy spectrum 
together with other previously shown [IJEJEl properties of it, imply that the bounded multibarrier 
potential is a singular system |H1 El EB EEU EE2 EE] • We conclude with a brief summary in Section 
6. 



2 The bounded one-dimensional multibarrier potential 

We consider a bounded one-dimensional system of N barriers assumed to have the same height V 
for all of them and to be uniformly arrayed between the point x = —^-^ and x = Thus, the 
total length of the system is L = a + b where a denotes the total width of all the N barriers (where 
the potential ^ / 0), and b is the total sum of all the intervals between neighbouring barriers 
(where V = 0). That is, in a system of ./V barriers the width of each one is and the interval 
between each two neighbours is j^t- Writing b as b = ac, where c is a real number, we can express 
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a and b in terms of L and c as [3] 



L 
1W 



T+c 



(i) 



We consider the passage of a plane wave through this system, which has the form 



A e 



1 k X 



+ 



B§e~ lkx where x < — Matching boundary conditions at the beginning and end of each barrier, 
we may construct a solution in terms of the transfer matrix ^ [21 E] ° n the J-th barrier. 
Thus, using the terminology in pQ, we obtain after the n-th barrier the following transfer matrices 
equation [I] 



p{n)p{n-X) _ _ _ p(2)p(l) 



A 
Bo 



V 



B 



(2) 



B211+I 

where A2 n +i and i?2n+i are the amplitudes of the transmitted and reflected parts [I] respectively 
of the wave function at the n-th barrier. Ao is the coefficient of the initial wave that approaches 
the potential barrier system, and Bo is the coefficient of the reflected wave from the first barrier. 
V is the total transfer matrix over the system of n barriers and is given for the E>V case in the 
limit of very large n by [I] 



Pe>v 



e lz (cos(j) + if 



sin((f>) - 



e lz (cos (j) — if sm ^ 
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The parameters f,d,z and (f> are expressed as jl] 



f = kb + aq^, d = aq^, z = k(a + b), = y/ f 2 - d 2 , 



where k, q, £ and r\ are [1] 

k = 



2m(E) 



2m{E - V) 
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It is shown 0] that V is given for the E<V case and in the same limit of very large n by 



(4) 



(5) 
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where the parameters /, d and (/> are now 



f = kb-^>, d= a 4, 

J 2 2 



/ 2 -d 2 , 



(7) 



/c is the same as for the E > V case and q = y 2m (^ g ) . i n the numerical part of this work we 

follow the convention in the literature that assigns to h and m the values of 1 and \ respectively. 
Defining, as in jSJ, the parameters k and r 



/sin((/>) 



cos(c/>) + V— 



d 2 sin 2 (c/>) 



r= 1 + 



(8) 



3 



one may obtain the eigenvalues of either the matrix of Eq (J2J) for the E>V case or of the matrix 
(jHJ) for E < V in the form 



Ai,2 = tcos(0 — «) ± \J t 2 cos 2 (4> — k) — 1 (9) 

The eigenvalues of the E>V case are obtained by substituting the <f> and / from Eq @ and those 
of the E<V case by substituting the corresponding quantities from Eq 0. 

We, now, show that, unlike the infinite Kronig-Penney system ^ [2J |3] that have no finite total 
length L and no finite ratio c, the energy spectrum of the bounded multibarrier potential depends 
critically upon L and c. We first take the real components of both sides of the first of Eqs (jHJ) and 
divide the numerator and denominator of its right hand side by cos(i^) 

cos(k) = 1 = (10) 

We note that the last equation is valid for both cases of E > V and E < V except that one should 
substitute the appropriate / and (f> from Eq © for the E>V case (as we do in next section) or / 
and 4> from Eq ((7J) for E < V as done in Section 4. 



3 The band-gap structure of the finite multibarrier system for the 

E>V case 



We, now, consider the E>V case and substitute in Eq (fl7)|) the appropriate /, d, (ft, k and q from 
Eqs (jlj) and © to obtain 



cos(k) 



1 + (1 + 4E(JWWW tan2(jB( ° + " V{a2 + ^ 

(ID 



i + (i + mi+c/m^-v) ) ^ 2 (l 2 (e - £)) 



The last result was obtained by expressing a and b in terms of L and c according to Eq lfl|). Squaring 
both sides of the last equation one obtains 



CO s 2 (k) = -5 rr — (12) 

1 + (1 + im+m+cFvT ) ^ 2 mE ~ &)) 



one must have corresponding values of L 2 (E — j^) that cause tan 2 (L 2 (E — j^pj)) on its right 
hand side to become very large. That is, for k = ± ( 2N ^ n ; (jy = 0,1,2,3,...) one must have 



From Eq l|12() we see that for all values of k that cause cos 2 (k) on its left hand side to vanish 

^) that cause tan 2 (L 2 (E - ^) 

[+. 
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I?{E — = ( 2jV + 1 ) 7r ; ^jy = o, 1, 2, 3, . . .). Note that since E, V and c are positive the expression 
L 2 (E — for the E > V case is positive. Thus, the appropriate values of the energies E that 
correspond to the remarked values of k are E = ^ 2 ^^P n + jxj. That is, from the last expression 
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for E and from the condition E > V we find that there are allowed energies only for those values of 
the total length L and the ratio c that satisfy the inequality 

( 2N + m + o >Vt N = Q1 ^ (13) 

That is, the energies that correspond to k = ± ( 2JV + 1 ) 7r ; (jy = 0, 1, 2, 3, . . .) depend upon the total 
length L of the system and the ratio c of its total interval to total width. For all other values of k 
we may take the reciprocals of both sides of Eq (|T2*|) and subtract 1 from the resulting expressions 
to obtain 

ton2(K) = (1 + 4 E(1 + c )(I(i + c )-K) )ton2 ' L2 ' E - TT^ (14) 

We note that exactly the same equations as those of (|12|) and (j!4jl are obtained also for the E < V 
case except that we use, as remarked, the / and 4> horn Eqs (J7J) with q = \J 2m %T ^- (where, for 

numerical purposes, we assign, as noted, h = 1, and m = ^). From Eq (|Fi|) we may find the energies 
which correspond to the values of k that cause tan 2 (ft) on its left hand side to vanish. These k's 
are k = ±Nir, (N = 0, 1, 2, 3 . . .) so using Eq ((Til) and the positiveness of I?(E — we have 
1?{E — yq^) = Nir from which one obtains E = ^ + yq^. From the last expression for E and 
from E>V we find that there are allowed energies that correspond to ft = ±Nir, (N = 0, 1, 2, 3 . . .) 
only for those values of L and c that satisfy 

^£^>V, iV = 0,l,2,... (15) 

That is, as for the case of the inequality (|13|) . the energies that correspond to k = ±Nir, (N = 
0, 1, 2, 3 . . .) depend upon the total length L of the system and the ratio c of its total interval to 
total width. Note that although 1?{E — = corresponds also to the former values of k that 
result in tan 2 («) = there is no value of the energy E, for the E > V case, that corresponds to 
L 2 (E — = 0. This is because from L 2 (E — ^r^) = one obtains E = which, for positive 
c, does not conform to E>V. 



Now, as may be seen from Eq (|14j) . the dependence of the allowed energies upon L, c and iV 
is not restricted only to these values of E that correspond to k = ± ( 2jV + 1 ) 7r or to k = ±Ntt, 
(N = 0, 1, 2, 3 . . .). Moreover, as will be shown in the following, the form of the allowed energy, as 
function of k, depends in a peculiar manner upon the ratio c. That is, if c is small, say < c < 3, 
then the expression that multiplies tan 2 (L 2 (£' — ^xj)) on the right hand side of Eq (|14|) is larger 
than unity and the solution of it for the allowed energy E as function of k can be obtained only 
numerically. In such case we find, as shown in the following, that these £^'s change periodically 
with k. If, on the other hand, c > 3 then the expression that multiply tan 2 (L 2 (£ - ^)) on the 
right hand side of Eq (|14j) is very close to unity. For example, for c = 3 and E ~ V this expression 
equals 1.020 and it tends fastly to unity for values of c > 3 or/and E>V. Thus, for c > 3 we may 
safely approximate Eq ((11)) as 

tan 2 ( K ) w tan 2 (L 2 (£ - —— )), (16) 

1 ~~ \~ c 

from which we obtain tan(L 2 (£^ — yx^)) = ±tan(ft). The last equation results in L 2 (E — yxj) = 
Ntt ± k, (N = 0, 1, 2, 3, . . .) from which we have 

E= N^±^ + V^ (# = 0,1,2,3,...) (17) 
L z 1 + c 
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As seen from Eq 1)17(1 not all values of Ntt ± k are suitable for the E>V case but only those that 
satisfy (Ntt ± k) > ^ , e c . For example, for V = L = 15 and c = 3 the permissible values of iW ± « 
that have corresponding energies are only those that satisfy (Ntt ± k) > 2531.25. 

For other values of c from the range < c < 3 the expression that multiplies tan 2 (L 2 (E — jx^)) on 
the right hand side of Eq ()14() is, as remarked, larger than unity for the E> V case and one must 
resort to numerical methods for solving Eq 1)14(1 for E as function of k. In this case the obtained 
energies are quasi-periodic where the character and form of the quasi-period depend upon the total 
length L of the system and the ratio c of its total interval to total width. We have emphasized the 
word quasiperiod since as one may realize from a close scrutiny of the appended figures the forms of 
the energy as functions of k are not exactly repeated over the k axis but differ in the small details. 
It is also shown for both cases of E > V and E < V that the quasi-periodic character of the energy 
as function of k emerges only for values of E that are close to the potential V (close from below for 
the E<V case and from above for E>V). For example, for E>V and V = 15 the corresponding 
energy is found to vary in almost periodic manner only between E = 15.5 and E = 19.5. This 
is because for higher values of E the expression that multiplies the tangent function on the right 
hand side of Eq ()14j) tends to unity even for very small values of c. In such case the situation is 
the same as that formerly found for large c (c > 3) where the appropriate equation to use is 1)16(1 
which leads to the linear expression 1)17)1 for the energy as function of k. For values of E that are 
close to the potential V one may obtain not only the remarked quasi-periodic forms of the energy 
as function of n but for very small c the allowed energies become constant as demonstrated by the 
horizontal lines in Figures 1-2. This is because for E ~ v and c ~ the right hand side of Eq 1)14)1 
becomes the indeterminate form §, so using L'hospital theorem ^H] we obtain for it the value of 1. 
In such case the value of E ~ V does not change with k. 

The quasi-periodic nature of the energy as function of n is demonstrated for small c in the form of 
half squares as seen from Figure 1 which shows nine different curves of E from Eq (|14[) as function 
of k. All the nine curves are drawn for the same values of V = 15 and L = 100 but for 9 different 
values of the ratio c = 0.2 • n, (n = 1,2,3, ...9). The different graphs in this figure and in Figures 
2-3 correspond to the values of c in a such a manner that the curves with the most pronounced half 
squares, which are shown at the upper part of these figures, fit the lower values of c. The curves 
which are characterized by a short and flat half squares, which are generally located at the lower 
parts of these figures, correspond to the higher values of c. The horizontal lines at the bottom of 
Figures 1 and 2 correspond, as remarked, to the constant value of the energy E ~ V. The remarked 
correspondence between the graphs of E and the appropriate c's, by which the graphs at the upper 
part of the figure fit the lower c's and those at the lower part of it fit the higher c's, applies also 
for the E<V case in Figures 4-6. 

As seen from Figure 1 any two neighbouring half squares are connected to each other by short 
horizontal lines where the energy is clearly discontinuous at the connecting points. That is, at 
these points the energy jumps in a discontinuous manner to form the squared parts where the 
smaller is c the larger are the corresponding jumps as may be seen at the higher part of the figure. 
For growing c the jumps become smaller and the corresponding half squares shorter as seen at 
the lower part of Figure 1. One may also realize that for the large value of L used in this figure 
the 9 different curves are shown mixed inside each other between the two values of E ~ 15.3 and 
E ~ 17.8. Plotting the energies E as functions of k for additional larger values of c, while keeping 
the same L as before, add more periodic forms (not shown) that are arrayed inside the formers 
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between the same two limits of Figure 1. 



We must note that taking into account the periodic character of the tangent function we can rewrite 
Eq (|14(1 in its most general form as 

*n»(„ ± Nn) = (1 + iE{1+c){ V E l + c) _ v) ) ^\L\E - JL-) ± m), (18) 

where N, N = 0, 1, 2, Since the parameter L occurs only under the tangent function at the right 

hand side of Eq 1)18(1 we may numerically solve it for E and find the same quasi-periodic dependence 
upon k for different values of L. For example, Figure 4, which was drawn for L = 30 may be obtained 
in almost the same form also for L = 110. Figures 1-6 of this work are obtained from the equivalent 
equation ([14)1 in which N, N = 0. Note that for very small L one may approximate the tangent 
function at the right hand side of Eq (|14jl by L 2 (E — and numerically solves the resulting 
equation for E as function of k. The spectrum obtained is shown in Figure 3 for the E>V case 
and in Figures 6-7 for E < V. 



As L decreases the different curves begin to be separated from each other and to occupy different 
sections of the ordinate axis. This is shown in Figures 2 and 3 which are drawn under exactly the 
same conditions and for the same values of V and c as in Figure 1 except that now L assumes 
the values of L = 5 for the curves of Figure 2 and L = 0.3 for those of Figure 3. The highest 
half squares of Figure 2 are drawn around the value of E = 19.6 and those of Figure 3 around 
E = 29. Note that the half squares of Figure 2 are each connected to the horizontal linear sections 
at more than two points (compare with Figure 1). This is especially demonstrated in the curves 
at the higher part of Figure 2 which correspond to the smaller values of c. Also, as in Figure 1 
the discontinuous jumps in the energy are larger for the curves that correspond to small c's at the 
upper part of Figure 2 and smaller for the larger c's at the lower part of it. The almost periodic 
curves of Figure 3 which are all drawn for the value of L = 0.3 are seen to be rather curved than 
squared and the points at which the energy jumps to form the pronounced parts of it are two for 
each quasi-period as in Figure 1. Also, the larger is the value of c at the lower part of Figure 3 the 
smaller become the corresponding jumps of E. 



4 The band-gap structure of the finite multibarrier system for the 

E <V case 



We, now, discuss the E < V case and begin from Eq (|12|) which is valid, as remarked, also for this 
case. We first find the allowed energies which correspond to the values of n that cause the left hand 
side of Eq (|12|) to vanish, that is, to k = ± ( 2Ar + 1 ) 7r ; (jy = 0, 1, 2, 3, . . .). As for the E> V case these 

energies are E = ± ^^jr^ + yq^; but now we consider also the E's which correspond to negative 
values of k so long as this does not cause the energies to become negative. That is, considering the 
negative values of k we obtain from the last inline equation for E and from E > that L, c and V 
should be related by the inequality ^ 2I ^P n < jj—, in which case L could not be very small. Thus, 
taking into account both positive and negative values of k and the condition J?<l'we see that 
there exist allowed energies only for those values of L and c which satisfy 

± l 2N + m + c > <Vt (JV=M!2) ... ) (19) 
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Figure 1: The nine graphs of this figure, which show the energy E from Eq (|14|) as function of k, 
are all drawn for the same values of V = 15 and L = 100 and for nine different values of the ratio 
c = 0.2 • n, (n = 1, 2, . . . 9). The curves that have pronounced half squares at the higher part of 
the figure fit the lower values of c and those that have short and flat half squares at the lower part 
belong to the higher c's. The horizontal line at the bottom corresponds to the constant value of the 
energy obtained at E ~ V. Each quasi-periodic half square is connected at its two sides by short 
horizontal linear sections. One may clearly see that the energy is discontinuous and undifferentiable 
at the vertical sides of each half square. 



S 




Figure 2: The nine graphs in this figure, which show the energy E from Eq (|14|) as function of k, 
are drawn under exactly the same conditions as those of Figure 1 except that now L = 5. One 
may see that the graphs are now separated from each other, especially those at the higher part 
of the figure which are clearly detached from the lower three graphs. Note that each half square, 
especially at the higher part of the figure, is discontinuous at several points and not only at its 
two sides (compare with Figure 1). The horizontal line at the bottom corresponds to the constant 
value E sa V. 
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Figure 3: In this figure, which is composed as in Figures 1-2 of nine graphs of the energy E as 
function of k, the total length L is further decreased for all the graphs to L = 0.3. The values of V 
and c remain as in Figures 1-2. Note that this additional decrease of L causes the almost periodic 
half squares to be narrower, more curved and arranged separately from each other compared to 
Figures 1-2. 
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In order to satisfy the last inequality for the case in which the minus sign at the left hand side 
is considered the parameters L, c and V must also be related, as remarked, by ^ 2I ^P n < 
Comparing Inequality 1)19|) with that of 1)13|) we see, as for the E>V case, that the allowed energies 
which correspond to k = ± ( 2jV + 1 ) 7r ; (_/y = 0,1,2,3,...) critically depend upon the values of the 
total length L and the ratio c. Note that whereas in the E>V case there are no allowed energies 
for the very large values of either L or c (or both, see Inequality <|13|) ). here, when considering the 
plus sign at the left hand side of (|19jl. we find no such allowed energies for the very small values of 
either L or c (or both). 

We, now, find the allowed energies that correspond to the values of k that cause the left hand side 
of Eq (fTl|) to vanish, that is, to k = ±Nir, (N = 0, 1, 2, 3 . . .). These energies are found, as for the 
E>V case, from L 2 (E — j^) = ±iW and are E = ±^f- + j^.- Using the former discussion at the 
beginning of this section we conclude that we may consider also negative k's so long as < ^ 
in which case L could not be very small. Thus, considering both positive and negative values of 
k and the E < V condition we find that there are allowed energies that correspond to k = ±Ntt, 
(N = 0, 1, 2, 3 . . .) only for L and c which satisfy 

±%^<V, (iV = 0,l,2,...) (20) 

In order to satisfy the last inequality for the case in which the minus sign at its left hand side is 
considered the parameters L, c and V must be related also by jj- < Comparing Inequality 
(|2U|) with that of (|15[). which refers to the E>V case, we find a result analogous to that formerly 
found from comparing the Inequalities (|13|) and (|19j) . That is, considering the plus sign at the left 
hand side of ()20|) . one may realize that for very small values of either L or c or both of them there 
are no allowed energies that correspond to k = ±Nir, (N = 0, 1,2, . . .). This is to be compared 
to the corresponding E > V case and the Inequality (|15|) from which one finds that there are no 
allowed energies for the very large values of either L or c or both of them. 

As realized from Eq l|14jl. which is valid also for the E < V case, the dependence of the allowed 
energies upon L and c is not restricted only to these values that correspond to 
to k = ±Nir, (N = 0,1,2,3...). Thus, one may find also the energies that correspond to other 
values of k but in this case we have first to exclude those values of E that cause the expression 
that multiply tan 2 (L 2 (.E — jj-^j) on the right hand side of Eq l|14|) to become negative. Otherwise, 
the right hand side of Eq (|14|) would be negative whereas its left hand side is positive. Thus, the 
allowed energies are those that satisfy 



V 



2 



AE(l + c)(E(l + c)-V) 



> "I (21) 



The last inequality results in 



V V 

E - E + — -7T > 0, (22) 

l + c 4(1 + c) 2 



from which one may conclude that the allowed energies of the bounded multibarrier system for the 
E<V case should be only those which satisfy the inequality 

V > E > (23) 
l + c 

Thus, the energies that correspond to k ^ anc | K _^ -j-jy^ should belong to the range 

V > E > which means that for very small values of c there are no allowed values of E at all. 
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For not very small c the energies are found by discussing, as in the E > V case, the two cases of 
large and small c. That is, for values of E that are close to V (but always E < V) the expression 
which multiplies tan 2 (L 2 (E — yq^;)) on the right hand side of Eq (|14|) tends to unity for c > 3. In 
such case we may use, as in the E>V case, the simpler Eq (j!6|l which results with the expression 
(|17[) for the allowed energies. But now since, as remarked, E satisfies V > E > there are 

corresponding energies only for < Nir ± k < (N = 0, 1, 2, . . .). Thus, using Eq (|17|). which 

is valid also for the E < V case, we find that the energies which correspond to c > 3 and which are 
close from below to V are allowed only for k, c, L and V that are related by 

('+■=>("' *«> < V (24) 
L z c 

For values of c from the range < c < 3 and for energies that are not close to V we use the full 
equation (|14[) which necessitates numerical methods for solving it for E as function of k. It is 
found, as in the E > V case, that the required energies are either quasi-periodic or constant. The 
constant energies are demonstrated by horizontal lines and appear when E assumes the value of 
E ps In this case the right hand side of Eq (|14jl assumes the indeterminate form R so using 
L'hospital theorem |16| it becomes unity for all values of k. The analogous cases discussed in the 
former section for E>V have resulted in finding constant horizontal lines for E ps v. 

Each of the following three Figures 4-6 shows 14 different graphs of the energy E as function of k 
for the same value of v = 15 and for c = 0.2 • n, n = 1, 2, 3, . . . 14. As for Figures 1-3 the 14 curves 
in each figure of the set 4-6 fit the 14 values of c in an inverted order. That is, the highest value of c 
fits the lowest curve in each figure and the second value of c from above corresponds to the second 
curve from below in each figure and so forth. The 14 graphs of Figure 4 are all drawn for the value 
of L = 30 and one may see how the energy, in the form of half squares, changes almost periodically 
with k. As for the E > V case the forms of the energy are not exactly periodic along the k axis. 
One may also realize that the most pronounced half squares are obtained for the smallest value 
assumed here for c (c = 0.2). Each half square is connected at its two sides to its neighbours 
by horizontal lines where, as for the E > V case, the energy jumps at the connecting points in a 
discontinuous manner. Note that the larger is c at the lower part of the figure the longer become 
the connecting horizontal lines between the half squares which become shorter and narrower. 

The 14 different graphs of the energy E as function of k, which are shown in Figure 5, are all drawn 
for the same values of V and c as those of Figure 4 but for L = 5. The squared sections are now less 
emphasized compared to those of Figure 4 and they become flattened for the larger values of c as 
seen in the curves at the lower part of the figure. Also, compared to Figure 4, each squared section 
is connected, especially at the higher part of Figure 5, to more than two points at which the energy 
jumps discontinuously. From Figures 4-5 and from the occurence of L 2 under the tangent function 
in Eq 1)14(1 one may realize that decreasing the value of L results in a corresponding decrease of the 
height and width of the quasi-periodic squared sections. This is clearly shown in Figure 6 which 
shows 14 graphs of E as function of k that are drawn under exactly the same conditions and the 
same values of V and c as those of Figures 4-5 except that now L assumes the rather small value 
of L = 0.8. The quasi-periodic parts, which are no longer half squares, have been considerably 
flattened to the degree of becoming almost horizontal lines as seen clearly in the curves at the lower 
part of the figure. 

Continuing to further decrease the parameter L results in widening the gaps in the energy spectrum, 
that is, in finding no corresponding allowed energies for entire ranges of the k axis. This is clearly 
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shown in Figure 7 which shows the form of E as function of k for the same V = 15 and the same 
values of c = 0.2 • n, n = 1, 2, 3, ... 14 as those of Figures 1-6 but for L = 0.278. As seen from the 
figure the allowed energies are represented by vertical lines between E = 5 and E = 15. These lines 
differ by their lengths which correspond to the different values of c. The thick lines signify that 
several lines which correspond to different values of c are identical and are seen tangent to each 
other. Note also that the vertical lines tend to be collected together (see, for example, the collection 
of the 14 lines at k, ~ 20) and this gathering is repeated almost periodically along the k axis where 
there exists no energy between any two neighbouring quasi-periods. That is, the gaps between the 
allowed energies are quasi-periodically repeated over the positive k axis and do not vanish for any 
value of it. This result denotes, as will be shown in the next section, that the bounded multibarrier 
potential is singular for this value of L. 

In Table 1 we show the allowed energy E and its dependence upon k, L, c and V for both cases of 
E>V and E<V. 

Table 1: The table shows the allowed energy E for specific values of k and for both cases of E > V 
and E < V. We may realize from the table that the allowed energies of the bounded multibarrier 
potential depend critically upon the total length L and the ratio c. Thus, one may conclude that 
for certain values of L and c there are no corresponding energies as discussed in details in Sections 
3 and 4. Also, there are no energies, especially, in the limits of very small L or c (or both) for the 
E < V case and in the limits of very large L or c (or both) for the E>V case. By N we mean all 
the natural numbers and zero (N = 0, 1, 2, 3, . . .). 



K, C 


E>V 


E<V 


For all k and small c 

_ (27V+1)tt 
— 2 
_ (27V+1)tt 
2 

k = Nir 

K = -NTT 

Ntt±k, k^±Ntt, k^± (2N ^ 1)7T , c>3, 
Nit±k, k^±Ntt, k + ± (2JV + 1V , 0<c<3, 


E constant for E ~ V 

j? (2N+l)n , V 
^ 2L^ + 1+c 
p (27V+1)tt , V 

& - 2L^ + 1+c 
_ Ntt , V 
~ ^ l+c 

p _ Ntt , V 

~ ' l+c 

h —[7- + 1+C 

(quasi-periodic, see text) 


E constant for E ~ 

j? (2N+l)rr , V 
^ 2L^ + 1+c 
tp (2N+1)tt , y 

^ - 2L^ + 1+c 
it. Wrr 1 y 

^ - + —c 

TP Ntt 1 y 

^ - ~TT + l+H 
y ^ tp Afrr+K 1 y ^ 

1+c < ^ T? + 1+c < K 

(quasi-periodic, see text) 



5 The singular character of the bounded multibarrier system 

We have seen in the former Sections 3-4 that the bounded multibarrier potential have gaps associ- 
ated with its energy spectrum which depend upon the total length L of the system and the ratio c 
of its total interval to total width. We also show that, unlike the band-gap structure of the infinite 
Kronig-Penney system in which the gaps disappear for large values of the energy, the remarked 
dependence here of the energy spectrum upon L and c implied that for certain values of them the 
gaps do not disappear. This has been shown for the E > V case in Section 3 and especially for 
E < V in Section 4 where we see that there are no energies that satisfy E < r±+£\ ( see Inequality 
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Figure 4: The 14 graphs of this figure show the energy E from Eq (|14|) as function of k for the same 
V = 15 and L = 30 but 14 different values of c = 0.2 ■ <n,(n = 1,2,... 14). The allowed energies 
are found in the range jt— < E < V. The graphs fit the values of c in an inverted order, that is, 
the higher c's correspond to the lower curves of the figure and vice versa. Note that for this value 
of L = 30 the almost periodic half squares are large and separated for the smaller c's at the higher 
part of the figure while they become narrow, flat and mixed among themselves for the larger c's at 
the lower part of it. 
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Figure 5: All the 14 graphs of this figure, which show the energy E from Eq (|14jl as function of 
k, are drawn for the same V = 15 and the same values of c = 0.2 • n, (n = 1,2,... 14) as those of 
Figure 4 but for the lower value of L = 5. As for Figure 4 the curves fit the values of c in descending 
order so that the lower values of c correspond to the curves at the higher part of the figure and 
the higher c's fit the curves at the lower part of it. Note that the quasi-periodic squared sections 
are less pronounced compared to those of Figure 4 and they become almost horizontal lines for the 
higher c's at the lower graphs. 



15 




4 - 



2 - 



2 4 6 8 10 12 14 16 18 



Figure 6: The 14 graphs of this figure, which show the energy E as function of n, are drawn under 
exactly the same conditions and the same V and c's as those of Figures 4-5 except that now L is 
decreased to the value of L = 0.8. Note that the curves become almost horizontal lines and this is 
further emphasized for the higher values of c in the lower graphs of the figure. Also, one may see 
that the discontinuity encountered for the larger values of L in Figures 4-5 is now less pronounced. 
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Figure 7: This figure shows the form of the energy E, for the E < V case as function of k, for 
the same V = 15 and the same values of c = 0.2 ■ n, n = 1, 2, 3, . . . 14 as those of Figures 4-6 
except that now L is further decreased to the value of L = 0.278. As seen, the allowed energies 
are represented by the vertical lines between E = 5 and E = 15 where the thick lines signify that 
several lines are identical and are drawn tangent to each other. One may realize that the gaps for 
this small L are repeated over the whole positive k axis and do not vanish for any value of it. 
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(J22J) and also Figure 7 from which we realize that there exist quasi-periodic gaps for small values 
of L). Now, it is accepted in the literature (see, for example, [H1E3]) that if the energy spectrum 
of any system have gaps which do not diminish for either large values of the energy E or large n 
then the relevant system is considered singular. We have shown in Sections 3-4 that not only the 
energy spectrum of the bounded multibarrier potential have gaps which remain constant along the 
K axis but that even for the parts of the epectrum which have bands there are points at which the 
energy, as function of k, is discontinuous. These points are shown, for both cases of E > V and 
E< V, to be associated with large L and small c as seen in the higher parts of Figures 1-2 and 4-5. 
Only at small L one may obtain energies that change continuously with k and do not jump in a 
discontinuous manner (see Figures 3 and 6). Thus, the existence of the points at which the energy 
changes discontinuously with k together with the remarked nonvanishing gaps demonstrate that 
the bounded multibarrier potential is singular. This may be corroborated by a previously obtained 
result [U |Hj which shows that the transmission probability of the bounded one-dimensional multi- 
barrier is unity. That is, the incoming wave function remains the same (preserved) after passing the 
dense array of potential barriers [HIS]- Thus, the three properties which are characterized in the 
literature jSJ |UJ 1101 lllj as signs of singular systems are all found also in the bounded multibarrier 
potential discussed here. These properties are : (1) nonvanishing gaps, (2) the ability to pass all 
the barriers of the system and (3) zero range potential for each barrier which arises here from the 
fact that the bounded spatial length contains very large number of barriers. Thus, we conclude 
that the dense array of potential barriers is also singular. 

A singular behaviour of the energy spectrum is demonstrated also in what is termed the singular 
point interaction 5 which is characterized by a periodic array of identical short-range perturbations 
[SJ 1101 lllj like our system. The allowed bands resulting from this 5 interaction are characterized 
10 by an increasing, as E grows, of the gap to band ratio. The difference between the periodic 
singular interaction 5 and our system is that the length of the former is infinite ^U] , whereas our 
system is bounded. 

We note that infinite spatial length is also what characterizes many examples of singular systems 
discussed in (Hj, so the physical interpretation of the singular continuity of these systems (see [9lll0|) 
is connected to their infinity. That is, as written in [S] "the particle will be transmitted at least once 
through each barrier and then it will return again to the initial point from which it has started". 
Thus, for these infinite systems the singular continuity is characterized by complete reflection that 
follows a former stage during which all the barriers are transmitted. This is so for the infinite 
systems in which no kind of transmission is possible at all, but since the bounded system admits 
transmission then the particle after passing all barriers finds itself outside the system at the other 
side of it and not at its initial side. Thus, what characterizes its singularity is complete transmission 
after the particle has been transmitted through each barrier. This was clearly shown jlj [3] , using 
transfer matrix methods P3E1E], for the bounded system discussed here. 

We, now, show this again by using the same analytical methods applied in ^U] for showing complete 
reflection of the 5 system. That is, taking into account the properties of the bounded multibarrier 
system and incorporating them into the formalism of ^U] we show that one obtains complete 
transmission instead of complete reflection. The single change we introduce into the formalism 
in ^U] is that of having to take into account that there exists only one kind of interaction in the 
bounded multibarrier potential compared to the multi-channels ITJ version of the 5 interaction 
discussed in [10 . That is, this 5 interaction is characterized as being composed from all possible 
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kinds of interactions (elastic, inelastic etc). Thus, the physical system used in ^U] is an infinite 
sequence of onionlike N channels scatterers, so that the reflection amplitude from any such scatterer 
is given in ^U] by 

-N 2 + 1 

R= N 2 + 2iN cot{kl) + 1' ^ 



where I is the length of each scatterer |10l and k = y ^p 5 - The total reflection character of 
this system is demonstrated |10j when the length I tends to zero and the number TV of channels 
(different interactions) in each scatterer becomes very much large. In these limits one may write 
jlUj cot(kl) ~ X and the unity term in Eq (|25|) may be discarded compared to N 2 . Thus, denoting, 
as in ^U], the product Nl by (3 the reflection amplitude from equation becomes JTCF R = 
— (1 + -fl) -1 - When, as in [Hlj, k — > oo (and still kl <C 1), the product (3k satisfies, due to the 
large N , (3k = kNl —* oo and the reflection amplitude tends to -1 so that the reflection probability 
R 2 becomes R 2 — > 1. Now, taking into account that for the single-channel multibarrier potential 
discussed here N is unity we trivially obtain, when substituting N = 1 in Eq R = 0. Thus, 
from the relation between the transmission T and the reflection R 1 T = 1 — R we learn that 
T = 1 as obtained from the transfer matrix method in We show in the Appendix that if the 

total length of the system is bounded then even if the interaction is of the multi-channel type, as 
in the 5 interaction, one obtains zero value for the reflection R provided the number of barriers n 
is much larger than the number N of channels. That is, if the total length L satisfies L < oo and 
if n » N then R — > which means that T — > 1. 

In summary, we conclude that if the total length of the system is finite then the singular character of 
the bounded dense array of either single or multi-channel scatterers is affected through the complete 
transmission it demonstrates. 



6 Concluding Remarks 



We have found the allowed energies of the bounded one-dimensional multibarrier potential for both 
cases of E > V and E < V. These allowed energies critically depend upon the values of the total 
length L of the system and the ratio c of its total interval to total width. Thus, it has been shown 
that there are values of L and c for which there are no allowed energies. For example, it has been 
demonstrated in Section 4 for the E < V case that the allowed energies are only those from the 
ranee V > E > -^- c . That is, for very small values of c there exist no energy for the bounded 
multibarrier system (see also the constant gaps in Figure 7). 

It has, also, been shown for both cases of E > V and E < V that for small c and large (or 
intermediate) values of L the energy E have the form of quasi-periodic half squares which are 
connected to each other by horizontal lines. The points of connection of these half squares, as may 
be seen from the figures, are points of discontinuity at which the energy can not be differentiated. 
For small values of L these squared parts become shorter, narrower and more curved (see Figures 3 
and 6) until they disappear for small enough L in which case the energy, as function of k, resembles 
the form of either horizontal (or vertical lines for the E < V case, see Figure 7). 

The remarked existence of gaps which do not disappear for certain values of L and c, even for large 
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k's, implied that the spectrum of the bounded multibarrier potential is singular. For example, it 
has been found for the E < V case and for very small L that the allowed energies have the form 
of vertical lines with nonvanishing gaps among them as seen in Figure 7. These gaps are quasi- 
periodically repeated along the k axis and do not vanish for any value of it. The constancy of gaps 
in the energy spectrum of other physical systems, like the 5 mentioned in Section 5, is interpreted 
in the literature |10| as resulting from the singularity of the involved systems. 

We note that the spectrum of the infinite Kronig-Penney periodic potential is absolutely continuous 
|13| . In such infinite periodic system one can not define, in contrast to bounded systems, any total 
length L or any ratio c of its total interval to total width. That is, the gaps in their energy 
spectrum, which do not depend on any such undefined L and c, can not be preserved, as in 
the bounded potential, merely by taking some limit or other values of these L and c. Thus, their 
spectrum is, as remarked, generally absolutely continuous |13| . There are, of course, the exceptional 
infinite periodic systems 111 j mentioned in Section 5 which are characterized by nonvanishing 
gaps but for other reasons not related at all to any L and c. The bounded multibarrier potential 
discussed here allows one not only to define its total length L and the related ratio of its total 
interval to total width but also to express the location of each barrier [IJ |SJ |S] and consequently 
the allowed energies in terms of L and c. Thus, there may be found, as actually shown in Sections 
3-4, specific values of L and c for which the gaps in the energy spectrum do not disappear even 
for large values of E or (and) large k. Moreover, for large (and intermediate) L and small c the 
energy E, as function of k, of the bounded multibarrier potential have been shown to have the 
form of quasi-periodic half squares which are jumpy and discontinuous at their vertical sides (see 
Figures 1-2 and 4-5). These discontinuities in the energy become small and insignificant only when 
L decreases as may be seen in Figure 3 for the E>V case and in Figure 6 for the E < V one. When 
L becomes very small the gaps in the energy spectrum becomes wider and are repeated along the 
positive k axis as seen in Figure 7. 

In summary, we see that the bounded one-dimensional multibarrier potential demonstrates through 
the discontinuities of E, the nonvanishing gaps in its energy spectrum, its almost zero-range po- 
tential and its complete transmission that it is a singular system. 
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A APPENDIX 

The transmission for the bounded one-dimensional multi-channelled 
multibarrier potential 



We, now, discuss the case of the bounded one-dimensional multi-channelled multibarrier potential 
system and show, using Eq (|25l and the relation R + T = 1, that the transmission tends in this 
case to unity The length of each scatterer (barrier) is now I = jj- where L, for the bounded system, 
signifies the total length of all the scatterers (barriers), and is some finite number, n is the number 
of scatterers in the bounded system and N is the number of channels in each scatterer where we 
assume n ^> N. We assume in the following that iV is large but small compared to n which, 
actually, satisfies n — > oo. Thus, denoting as before the product of the length of each scatterer and 
N by j3 (so as to have now IN = = j3) and using the former approximations of cot(kl) ~ h 
and N 2 + 1 « N 2 , we may write Eq ((2*5]) as 

-N 2 + l -N 2 + l , , 

N 2 + 2iN cot(^) + 1 ~ N 2 + + 1 ~ 

-N 2 2i x 

N 2 (l + ^)~~ [1 + W 

Now, when k — > oo (and still kl = — <C 1 as for the analogous case of Eq (|25|)) we have now, 
due to n S> N, k/3 = MdL <g \^ an( ^ ^j^g reflection amplitude goes to zero in contrast to the result 



obtained from Eq "25") . Thus, from the relation between the reflection and transmission amplitudes 
T = 1 — R, we learn that T goes to 1 as required. That is, we see that if the length of the sequence 
of scatterers is finite and if it is dense enough so that the number of scatterers n is very large then 
the transmission coefficient goes to 1 even for large number of channels N in each scatterer. 
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